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Abstract 

Following the recently accumulated information on the vibrational proper- 
ties of the Cu(OOl) surface acquired through single- and multiphonon He atom 
scattering experiments and the concomitant theoretical investigations, we have 
reexamined the properties of the Debye- Waller factor (DWF) characteristic 
of the He— >Cu(001) collisions using the recently developed fully quantal and 
three-dimensional model of inelastic He atom scattering from surfaces. We have 
focused our attention on the role which the various He-surface model poten- 
tials with their characteristic interaction parameters (range of the interaction, 
momentum and energy transfer cut-offs etc.) employed in the interpretation 
of the scattering data may play in determining the magnitude of the DWF. 
By combining the He-Cu(OOl) potential whose repulsive and attractive compo- 
nents are both allowed to vibrate with the substrate phonon density of states 
encompassing anharmonic effects, we obtain the values of the DWF which agree 
nicely with the experimental data without invoking additional fitting parame- 
ters. On the other hand, by taking the phonon momentum transfer cut-off Q c 
as an adjustable parameter, as has been frequently exploited in the literature, 
all the considered potentials can produce agreement with experiments by vary- 
ing Q c . The magnitudes of such best fit Q c values are compared with those 
available in the literature and their physical significance is discussed. 



"Contact author. Permanent address: Institute of Physics of the University, P.O. Box 304, 10001 
Zagreb, Croatia. E-mail: branko@ifs.hr 



1. Introduction 



A number of experimental studies of the vibronic properties of single crystal sur- 
faces carried out in the past decade by utilizing thermal energy He atom scattering 
(HAS) have emphasized the importance of the microscopic properties of the interac- 
tions between He beam atoms and the surfaces investigated. The interpretation of 
both the single phonon HAS data and the multiphonon scattering spectra required 
a relatively detailed specification of the form of the interaction potentials and the 
corresponding matrix elements. The latter proved to be one of the key quantities in 
establishing a meaningful comparison between experiment and theory [|T], ^ |3], |], |5|, |J . 

In a number of investigations of the vibrational properties of solid surfaces a suc- 
cessful theoretical interpretation of the relative single-phonon intensities was achieved 
by using adatom-substrate pairwise potentials and the distorted wave Born approxi- 
mation (DWBA) for description of the projectile motion in the static component of 
this potential [I], [3j. The studies of Cu(OOl) surfaces by HAS emerged as a partic- 
ularly illustrative example in the context of the atom-surface interaction potentials. 
To interpret the He— >Cu(001) time of flight (TOF) spectra certain refinements of 
the theory based on the pairwise potentials @ have been proposed either through 
the concept of a pseudocharge model |Q, H or the introduction of anisotropic He- 
substrate atom pair interactions 0. In both cases the modified matrix elements 
were typified by some cut-off parameters whose variation could significantly affect 
the one-phonon intensities. However, thus calculated one-phonon scattering spectra 
could reveal correctly only the relative intensities of distinct phonon modes. 

An additional test of the accuracy of some aspects of the scattering potentials 
may be provided by the multiphonon spectra whose dependence on the details of the 
interaction is more complex. Namely, in the multiphonon collision regime already 
small changes in the projectile-phonon coupling may give rise to a cumulative effect 
in the scattering intensity. Surprisingly enough, the intensities of the multiphonon 
scattering spectra of the He— >Cu(001) collision system || |[ were relatively success- 
fully reproduced theoretically || [TI| by using the earlier expressions for He-surface 
potentials and two somewhat different approaches to multiphonon HAS. How- 
ever, no consensus on the values of the characteristic potential parameters has been 
achieved in these two studies as in one of them || different sets of the interaction 
parameters have been introduced for single and multiphonon scattering regimes. 

Another basic quantity characteristic of the atom-surface scattering spectra which 
is sensitive to the features of the atom-surface interaction potentials is the so called 
Debye- Waller factor (DWF). In the surface scattering experiments it gives a measure 
of the intensity Io of the elastically scattered specular beam relative to the incoming 
beam intensity Jj„ and in this respect it differs from the notion introduced in neutron 
scattering from crystals fTI , p| . So defined DWF, commonly written in the form 



lol hn — exp(— 2W) where 2W is the corresponding Debye- Waller exponent, becomes 
essential in the multiphonon scattering regime because it provides the normalization 



proper of the scattering spectrum [U| [14], |15(] which should obey the unitarity prin- 



ciple (optical theorem). Also, in contrast to the measurements of the single phonon 
scattering spectra, the measurements of the DWF provide values which in a sense 
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are "absolute", i.e. unaffected by the time of flight (TOF) technique, angular and 
kinematic restrictions, etc. 

A rough estimate of the magnitude of the Debye- Waller attenuation in atom- 
surface scattering was given long ago by Weare [|ljj and later rederived by Levi and 



Suhl [11]. They arrived at an approximate expression for the DW exponent as a 



function of the substrate temperature T s in the form: 

Um 2W (T,) = 3 W fF) = 24 M "* E ' ms2e > I T - 



T s >0d " M crys k B QD V0D/ M crys k B QD \<S>d. 

Here Qjj is the surface Debye temperature of the substrate, Ak z is the change of the 
projectile momentum normal to the surface, Ei and $i are the incoming energy and 
angle of scattering of the projectile, respectively, M crys is the mass of the crystal atoms 
and ks is the Boltzman constant. Ei in this expression is sometimes also corrected 
for the surface potential well depth D (Beeby's correction [17|) in which case (hAk z ) 2 
is replaced by (hAk z ) 2 + 8Mff e D. However, the form of the DW exponent (Jj) can be 



justified only in the regime of impulsive scattering [IT], [12| and therefore its validity 
is of limited range. In particular, for incident energies typical of thermal energy 
He atom scattering from surface phonons and soft projectile-surface interactions the 



approximation of impulsive scattering has been shown to become unreliable [|J| for 
making quantitative comparisons with the experimental data. 

The measurements of the DWF in HAS from metal surfaces have been systemati- 
cally performed for Ag(l 11) [0, Pt(lll) @, Cu(OOl) ^Dj, Cu(llO) |HJ and Ni(115) 



[ p!q[ and theoretical interpretations given within two different approaches developed 
to treat multiphonon excitations is atom-surface collisions. The calculations of the 



DWF in HAS from Ag(lll) |23|] and Pt(lll) [19 were based on a three-dimensional 
scattering formalism developed in Ref. fl24| . They reproduced successfully the ex- 
perimentally observed magnitude and the linear ^-dependence of the Debye- Waller 



exponent 2W up to T s = 700K with the "vibrating soft atom model potential" |£3 
The experimental data on the DWF available for Cu(OOl) and Cu(llO) surfaces were 
interpreted by carrying out perturbation expansion of the scattering matrix in a dis- 



torted wave basis Although these calculations were essentially one-dimensional 



as regards the scattered particle dynamics, their real merit lies in finding that for the 
atom-phonon coupling to all orders in the lattice displacements the repeated single 
phonon exchange processes give much larger contribution to the DWF than the si- 



multaneous many-phonon exchange of the same multiplicity ||26|| . Upon introducing 
some adjustable parameters the DWF was calculated by assuming only the repulsive 
component of the total potential to vibrate and retaining in the scattering matrix the 
lowest order dominant contributions in powers of T s . Such a truncation of the series 
for the DWF, which violates the unitarity of the scattering spectrum, gives rise to 
an artifact in the curvature of the DWF versus T s on a semilogarithmic plot. 

Recently developed fully quantum, three dimensional (as regards the particle 
dynamics) and unitary multiphonon He atom scattering formalism |T5] embodies 
also a general expression for the DWF as an essential ingredient. The application of 
this formalism to multiphonon He— >Cu(001) collisions produced a nice agreement of 
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theoretical predictions |TIJ with the experimental multiphonon spectra || by using 
the same interaction parameters as in the one-phonon theory, i.e. different from 
those used by the authors of Ref. || to fit the experimental data. Hence, the state 
of affairs regarding the forms of the interaction potentials in HAS from Cu(OOl) 
surfaces seems to be critical as different expressions and parameters have been used 
to optimally describe the same physical situation. 

The selection of an adequate model atom-surface potential to describe inelastic 
processes in He— >Cu(001) collisions can be also appropriately tested in the evaluation 
of an "absolute" quantity such as the complete DWF provided the latter is calculated 
within a genuine three dimensional scattering model and by taking into account all 
dominant multiphonon contributions in powers of T s . This task is carried out in the 
present work by making use of the quoted novel multiphonon scattering formalism [|I3] 
[14], |l5l, [T0|. In the following sections we focus our attention specifically on the problem 
of which effects the various model potentials employed in recent interpretations of 
the one- JJ, |3], f|, |5], an d multiphonon HAS spectra from Cu(OOl) |9|, [TIJ may have 
on the magnitude of the DWF and its variation with T s . 



2. Atom-surface scattering potentials and the DWF 



In a series of earlier publications JT3|, [TJ], [15], |10[ we have shown that in in the case 
of linear projectile-phonon coupling, which gives by far the largest contribution to 
the multiphonon processes ||26|| , one can derive a closed form unitary expression for 



the scattering spectrum valid both in the one- and multiphonon scattering regimes. 
The point of departure in this approach is a distorted wave basis set of projectile 
wave functions which exactly describe elastic reflections from a flat surface. Inelastic 
processes are then treated as a perturbation of the distorted waves brought about by 
the vibrations of the surface. In the weak coupling limit the thus calculated inelastic 
scattering spectrum reduces, up to a kinematic Jacobian-like factor, to the standard 
inelastic reflection coefficient calculated in the DWBA [14], |HJ. On the other hand, 
for strong coupling the scattering spectrum is expressed in terms of uncorrelated 
and correlated multiphonon scattering processes in which only one phonon can be 
exchanged at a certain instant, and the normalization of the entire spectrum is given 



by the DWF ||15| . In the collision regimes typical of HAS experiments the uncorrelated 
multiphonon processes are dominant over the correlated ones and in this limit the 
scattering spectrum acquires a form of the exponentiated Born approximation (EBA) 



expression |T3], [14|, [15| . The corresponding DWF to all orders in the coupling constant 



reads 15 



e-™ EB A =eM-ERf^ B % (2) 

where Rf^ VBA is the temperature dependent one-phonon inelastic reflection coefficient 
|IJ, |2|, calculated in the DWBA. R^ BA is quadratic in the projectile-phonon 
coupling and describes the transition of the collision system from the initial state | i), 
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characterized by the particle distorted wave quantum numbers ki = (Ki, k iz ) and the 
initial phonon distribution, to a final state | /) in which the particle quantum numbers 
are kf = (Kf, kf z ) and the final phonon distribution differs from the initial one by 
one phonon quantum. Here frK denotes the lateral particle momentum and hk z the 
perpendicular particle momentum outside the range of the static distorting potential 
U(z) of the planar surface. Explicitly, one has [O, [UJ: 



2W eba = £ J, y£.Q .(+) | 2 [n ph {hu QJ ) + 1]+ I Vgfi^C-) | 2 n ph {hu Q ,)] (3) 

where Q, j and uqj denote the phonon wavevector parallel to the surface, the branch 
index and frequency, respectively, and n p h(^Q,j) is the Bose-Einstein distribution of 
phonons in thermal equilibrium at T s . The on-the-energy-shell one phonon absorption 



and emission matrix elements are given by 0, [14| 



VEfijCT) = ^vK£'f5 K , tK±Q 5(EK>, k > - E K , kz T huqj) (4) 

where E- K kz denotes the particle energy. | V^J3j(=f) 1 2 re P res ent first order DWBA 
probability for a state-to-state transition | K, k z ) — >| K', k' z ) of the particle. Their 
seemingly divergent contribution to (|3|) disappears upon conversion of the ^-functions 
to Kronecker symbols and summation over Q and k z . The matrix elements V kz ' k t j 
of the projectile-phonon interaction V(r) are taken with respect to the distorted 
waves corresponding to the projectile motion in U(z). A detailed derivation of these 
formulae and their application to HAS problems has been given in ]T5|, W, 



It is evident from Eqs. (|j) and @ that the DWF in atom-surface scattering 
will be sensitive to the form and variations of the dynamical atom-surface potential 
V(r) through its matrix elements V^'^j. These potentials are neither known a 
priori nor readily available in analytical form but have to be determined from inde- 
pendent calculations using various approximate schemes, often yielding them only 
numerically. On the other hand, in the three dimensional numerical calculations 
it is often convenient to have analytical representations of both the static and dy- 
namic He atom-surface interactions as this greatly simplifies the computing. This 
has stimulated the development of the various approximate analytical expressions for 
the potentials and their matrix elements as functions of the characteristic interaction 
parameters such as the strength and the range of the potential, typical cut-offs etc. 

In the majority of theoretical interpretations of inelastic one-phonon HAS exper- 
iments on smooth metal surfaces the static projectile-surface interaction is conve- 
niently represented by a sum of pair potentials v (r — ri) acting between the He atom 
at r and substrate atoms at rj: 

U(p,z) = J2<r-n), (5) 
l 

and then averaged over the surface [0] to yield the static projectile-surface interaction 
potential U(z). 
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The matrix elements of the dynamical interaction in the case of linear atom- 
phonon coupling acquire a simple form [|, ||]: 



V S'f = E u K (Q, j)F K (K - K', k z> , A; 2 )5 K -K', Q+ G, (6) 

G,K 

where the sum ranges over all reciprocal lattice vectors G of the surface mesh and the 
positions r K of the crystal atoms in the surface unit cell. u K (Q, j) is the quantized 
displacement of the crystal atoms corresponding to the phonon mode characterized 
by the quantum numbers (Q, j) and the polarization vector e re (Q, j). The matrix 
element of the force F K (Q, k z /,k z ) exerted on the projectile atom is expressed as 

F K (Q,k' z ,k z ) = e^( Xk ' z | -(iQ,^H4Q,z) I XO, (7) 

where | \k z ) is the distorted wave describing the particle motion normal to the surface 
and v V ib(Q,, z) is a two-dimensional Fourier transform of the vibrating part tw>( r ) of 
the pair interaction v (r). As yet, there is no general consensus on which part of the 
total pair potential contributes to tw>( r ) an d several models have been proposed in 
the literature. 

Another important parameter characteristic of the various expressions for v (Q, z), 
and thereby also of v V ib(Q,z), is the cut-off wavevector (or wavevectors) Q c which 
gives an approximate upper bound of the lateral momentum transfer in the one- 
phonon exchange processes (Hoinkes-Armand effect J30|). In the case of exponential 
potentials of range 1/(5 the value of Q c is approximately given by |7], |27j : 



(8) 

z t 

where z t is the average value of the He atom turning point in the surface potential 
U(p, z). Numerical evaluation of the matrix elements (0) avoids the explicit introduc- 
tion of Q c || but the physical effect of the cut off in the space of lateral momentum 
transfer persists. 

Expressions ©-© provide a framework for a fully three dimensional calculation 
of the DWF and thereby enable a test of the adequacy of the various expressions for 
the interaction potentials employed in the derivation of the matrix elements (|7|). 

The characteristics of the static He atom-Cu surface potentials have been ex- 



tensively discussed in the literature [31, 0, [33| 34, |35|, B6fl. Quite generally, this 



interaction is repulsive at short distances and exhibits an attractive van der Waals 
(vdW) tail at large distances, with a shallow potential well (6-7 meV) whose mini- 
mum is located around 7 bohrs away from the last crystal plane. The dynamic He 
atom-surface potential in the case of linear atom-phonon coupling is obtained as a 
gradient of the vibrating part of the pair potential, as implicit in Eq. (0). However, 
there is no unanimous agreement as to which part of the total potential this gradi- 
ent should involve, which is equivalent to the question whether only the repulsive 
or both the repulsive and attractive components of the total potential vibrate. We 
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shall calculate the force matrix elements (|7|) by using several different forms of U(z) 
and v V ib(Q, z) employed recently in the interpretation of both single and multiphonon 
scattering data on He— >Cu(001) collisions and then test their relevance by making a 
comparison with the experimental values of the DWF. 

3. The effect of the potentials on the DWF 

In our assessment of the effects of the various interaction potentials on the DWF 
pertinent to He— >Cu(001) collisions we shall investigate all three possibilities of dif- 
ferent forms of the static and dynamic interactions. Our point of departure will 
be the earlier calculated static He-Cu(OOl) potential |35j which was also in a good 
agreement with empirically determined potentials and experimental fits |j3i| . By re- 
quiring that the surface averaged sum of model pair interactions produces as close 
as possible this potential we may in principle determine the characteristic potential 
parameters corresponding to U (z) and v v n,(Q, z) which both derive from v(r). Then, 
by assuming various forms of iw>(Q, z) (originating either from only the repulsive or 
both the repulsive and attractive components of v(r)) we may be able to select the 
physically relevant U(z) and v v n,(Q,z) after comparing the calculated and measured 
values of the DWF. To pursue this goal but also remain in correspondence with earlier 
treatments of single- and multiphonon spectra in obtaining the wavefunctions needed 
for the calculation of the matrix elements (0) we shall fit the z dependence of the 
potential of Ref. |B5| to the following analytical expressions: 

(i) Exponentially repulsive potential 

U(z) = U exp (z) = U Q e- pz (9) 

by requiring that the value and the slope of the two potentials at the He-atom turning 
points be the same. This implies the pair potential in the form u(r) = Voe~ /3r and 
renders Uq and j3 as functions of the incoming energy Ei and angle of incidence 9i 
of the projectile. This procedure also enables us to estimate the values of Q c using 
Eq. (§) because the latter can be most easily justified in the case of exponentially 
repulsive surface potentials |7j. Here the whole potential is assumed to vibrate, i.e. 
v vi b(p, z) = v(p, z), and the calculation outlined in Ref. [7J leads to: 

«S?(Q, *) = U e-^ 2 ^ (10) 
with Q c given by @. The nontrivial component of the interaction matrix elements 
in (0), i.e. {xt^ P I — (*Q> ~§z) v vib (Qi z ) I Xk P )i are then obtainable in analytical form 

MM- 

(ii) The Morse potential 

U(z) = U M (z) = U™ p (z) + U^(z) = D(e- 2a{z - zo) - 2e-^ z ~ z ^) (11) 

where the well depth D, position zq of its minimum and the range a are determined 
by requiring that for given Ei and 0, the potentials (|TT|) and the one computed in 
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Ref. |35| have the same values at the minimum and at the classical turning point. 
An alternative requirement that the values of the minima and the derivatives at the 
turning point coincide leads to practically the same fitted Morse potential. The pair 
potential which leads to (|TTD is of the form v(r) = v (e~ 2ar — e~ ar ) and the different 
range of the repulsive and attractive components gives rise to different Q-dependence 
of v(Q, z). Here we assume that both components of the pair potential v(r) vibrate 
which leads to 

t# 6 (Q, z) = D( e -M*-*>) e -Q?/x% _ 2e-^-*°) e -Q 2 /Q?) (12) 

with Q c = J2a/zt- Here the factor of 2 in the second term in the bracket on the 

RHS of (|l|) and the different cutoffs in the Q-space (Q c versus Q c / \/2) arise as a 
consequence of the different range of the two components of v (r). Again in this case 
the corresponding matrix elements (xw I ~~ (*Q? ~jk) v vib(Qi z ) I Xkl) are available in 
analytical form [[39. |38] . 

(iii) The static potential U(z) is given by the Morse potential as in (ii) but only the 
repulsive component of v(r) is allowed to vibrate. This yields: 

<f (Q, Z) = De -Mz-z ) e -QV2Q^ (13) 

and the corresponding matrix elements are also expressible in analytical form [[f(], |38|j . 

The interaction potentials (|10D , (|i~2"|) and Cp~3f) and the corresponding matrix ele- 
ments should be representative enough to span the range of possible but physically 
different dynamical regimes of the projectile- surface encounter which the DWF can 
be sensitive to. The calculations were performed by using these potentials and the 
Debye model of surface phonons in Cu parametrized in terms of the surface Debye 
temperature 0d = 290 K which is a mean over the surface projected modes and 
directions in the surface Brillouin zone. This effective value is a little higher than the 
value of 267 K reported in || for the Rayleigh wave and the longitudinal resonance 
in (100) and (110) directions of the Cu(001) surface, and the value of 280 K found 
in |J as appropriate to the regime of multiphonon scattering from the same surface. 
The difference between the present and the other two temperatures arises because 
the latter were determined from fitting the strength of the projectile-surface coupling 
using a potential in which the effect of the attractive well was neglected. 

Lattice vibrations of copper surfaces exhibit anharmonicity even at moderately 
high temperatures and in order to take this effect into account we have additionally 
corrected the surface phonon density of states for anharmonic effects as discussed in 



Refs. [O, 42]. The actual effect of anharmonicity on surface phonon dispersion within 



the two dimensional Brillouin zone pertinent to the Cu(001) surface was estimated 



from the data presented in Ref. [ 43[ . 



In Fig. 1 we show a comparison of the results of calculations of the DWF for the 
three potentials described above, for £",=63 meV and 9i = 39° for which the shadow- 
ing effect in scattering from defects should not be very important. In addition to this 
we have also shown for illustration the DWF calculated by using Eq. (|l|) with the 
same 0£> as in other expressions. As is seen from Fig. 1, the Morse potential with 
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both components vibrating gives an excellent agreement with the experimental data. 
The Morse potential with only the repulsive component vibrating and the adjusted 
exponentially repulsive potential do not produce results which would adequately de- 
scribe the experiments for the given set of collision parameters, and neither does 
expression ([J). 

Fig. 2 is analogous to Fig. 1 in that it displays the same comparisons, but here for 
Ei = 21 meV and Q{ = 31.8°. The same general trends as observed in Fig. 1 persist 
for this set of collision parameters as well although the agreement is not as good as 
at Ei = 63 meV. Presumably this is due to the presence of diffuse elastic scattering 
which is more pronounced at lower collision energies but hasn't been accounted for 
by our model. A relatively good agreement between the measured DWF values and 
those calculated from Eq. ([]]) is a mere coincidence which doesn't occur at other 
collision energies and scattering angles (c.f. Fig. 1). 

A poor description of the DWF in terms of the adjusted exponential potential 
(designated by fllOD above) in Figs. 1 and 2 signifies the importance of the presence 
of the potential well, and this becomes more apparent as the normal component of the 
projectile incoming energy is lowered. On the other hand, the large difference between 
the results for the same Morse potential but with different vibrating components 
arises from smaller interaction matrix elements in the case of the total vibrating 
potential. This is mainly because the derivative of the complete potential changes 
sign at the well bottom giving rise to cancellation effects in the scattering amplitude 
(0). Another point worth emphasizing is that the lateral momentum cut-off function 
associated with the attractive component of the total vibrating potential attenuates 
much faster than the one associated with the repulsive vibrating component, giving 
correspondingly smaller contribution to (|7|). This simply reflects the fact that the 
potentials of longer range give rise to smaller inelastic contributions for given energy 
and momentum transfer. The same trends regarding the three model potentials are 
also retrieved for other scattering angles at He incoming energies of 63 and 21 meV. 
Hence, in the scattering regime characterized by the present collision parameters all 
these findings render, out of the potentials fllUp, ( |T2"D and (JIB|), the total vibrating 
potential designated by (|T2D as the best model potential underlying the physics of 
the DWF in the present collision system. 

In Figs. 3 and 4 we present further comparisons of the available experimental 
values of the DWF |20j with the ones calculated using the potential ([12]). Here we 
want to reiterate that in these calculations we haven't introduced any free parameters 
which could be adjusted so as to obtain a good fit to the experimental data but that 
all the parameters typical of the potential and the phonon density of states have been 
either predetermined or obtained within the present model. 

Fig. 3 displays the results of the calculations for incident energy Ei = 63 meV 



and the angles of incidence used in the experiments |20[ . The agreement between the 
calculated and experimental results is very good in the whole range of experimental 
8i in which also the earlier empirical fits were successful ||20|| . An exception occurs 
only at 9i = 71.6° for which due to the low normal component of the incident energy 
the contribution to scattering by surface defects (steps, kinks, adatoms, vacancies 
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etc.) may become important. In this situation two effects come into play. First, 
the diffuse scattering by static defects reduces the DWF at incident angles nearer 
to grazing because of the shadowing effect and, second, the enhanced anharmonicity 
of vibrations associated with defects would further suppress the magnitude of such 



The general trend of reduction of the DWF with the angle of incidence closer to 
the normal signifies the strongest coupling of the He atom to perpendicular surface 
vibrations. A deviation of the experimental results from this general behaviour for 
9i = 19° and 400 K< T s <800 K is probably an artifact connected with the early 



Fig. 4 displays the experimental and theoretical DWF values for incident energy 
Ei = 21 meV. Although the trends observed here are the same as in Fig. 3, due to the 
lower incoming energy the deviations between the two sets of data becomes apparent 
even at lower incoming angle. For $i = 73.5° the corresponding Ef is already so 
small (1.7 meV) that the scattering from defects may give rise to contributions to the 
Debye- Waller exponent 2W which are larger than that induced by phonons. As our 
model doesn't encompass this type of effects, the difference between the experimental 
and calculated values of the DWF for such collision parameters is not surprising. 

A crucial element in obtaining a good agreement between the calculated and ex- 
perimental values of the DWF without invoking the fitting parameters was a realistic 
form of the He-surface interaction potential with both components allowed to vi- 
brate, and the fact that the attractive component is much more sharply cut off in 
the Q-space due to its longer range. Thereby the attractive component gives only 
a correction to the one-phonon matrix elements (|?P whose magnitude is dominantly 
determined by the repulsive component of the interaction. 

Quite generally, the matrix elements ([/]) are so sensitive to even small variations 
of Q c that the latter may be taken as a fit parameter which could be adjusted so as to 
produce a good agreement between the calculated and measured values of the DWF 
for any of the potentials discussed under (i)-(iii) above. Fig. 5 shows a comparison 
of the magnitudes of such best fit Q c values together with the values available in the 
literature. However, a physical justification of such best fit Q c 's may not be easy in 
all the cases considered because some of them considerably deviate from the values 



predicted by the expression J2a/z t here found to provide a consistent description of 
the experimental results for the DWF in He— ^Cu(OOl) collisions. 

4. Conclusions 

In this work we have studied the effect of the characteristics of the various 
model interaction potentials on the magnitude of the Debye- Waller factor (DWF) in 
He— >Cu(001) scattering. For this collision system the experimental data on the DWF 
[p0|1 , single phonon ||, |, ||] and multiphonon spectra 0, || are available for a wide 
range of collision parameters (incident energy and angle) and surface temperature 
which all facilitates a comparison between the theoretical predictions and experimen- 



DWF. 



measurements . 
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tal results. In our theoretical description we have employed the earlier developed 
realistic, fully three-dimensional quantum scattering model to calculate the DWF to 
all orders in the projectile-phonon coupling [ 15| , here assumed to be linear in phonon 
displacements. The calculations were carried out within the so called exponentiated 
Born approximation (EBA) which takes into account uncorrelated multiple phonon 
exchange processes which have been shown earlier to give a dominant contribution 
to the scattering spectra in the collision regime studied [13|, |14], By considering 
several types of He-surface interaction potentials we have shown that the model po- 
tential which correctly reproduces the gross features of the earlier calculated static 
He-Cu(OOl) interaction and whose both components, repulsive and attractive, 
are assumed to vibrate, produces results which nicely agree with the experimental 
data without invoking any fitting parameters. Important elements in obtaining this 
agreement were the difference in the cut off wavevector Q c characterizing the longer 
range attractive and shorter range repulsive components of the vibrating potential 
and the variation of the phonon dispersion with temperature due to the anharmonic- 
ity of surface vibrations. Other model potentials which have also been frequently 
used in the interpretation of the DWF and the single phonon data, like the vibrat- 
ing exponentially repulsive potential, fail to reproduce the experimental values of 
the DWF in the He-^Cu(OOl) collision system. This is mainly because they do not 
take into account the attractive component of either the static or dynamic potentials 
which become important at lower incoming energies. On the other hand, at very low 
normal incoming energies the scattering from defects starts to affect the magnitude 
of the DWF. In this limit the agreement between the measured values and those 
predicted by the present model which doesn't account for scattering by defects, is 
lost. 

We have further shown that it is in principle possible to reproduce the measured 
DWF data with all the model potentials considered by assuming the characteristic 
cut off wavevector as an adjustable parameter for each set of collision parameters. 
Thereby one can determine the best fit Q c values as function of the normal component 
of the projectile incident energy for the studied model potentials and correlate them 
with the values given in the literature 0, |, [| [|, The best fit Q c values for the 



potential with both vibrating components were found to be practically identical to 
those of our no-fit-calculations, thereby confirming the consistence of this potential 
and the validity of the results obtained by using it. They are also found to be very 
close to the majority of the values cited in the literature in connection with the 
theoretical interpretation of both single- and multiphonon scattering spectra which 
gives additional credibility to these and the present calculations. 
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List of Figures 



Fig. 1. Debye- Waller factor as function of the substrate temperature for He incom- 
ing energy E{ = 63 meV and incident angle 9i = 39° calculated with three different 
interaction potentials. Full squares denote experimental datapoints [^0) and the long 
dashed, full and dashed-dotted lines denote the values calculated with the adjusted 
exponentially repulsive potential ([H])) total vibrating Morse potential ([Tj|) and re- 
pulsive vibrating Morse potential flT3p, respectively. The short dashed line denotes 
the DWF calculated from Eq. (§. 

Fig. 2. Same as in Fig. 1 but for E, L = 21 meV and 0* = 31.8°. 



Fig. 3. Comparison of calculated (full lines) and experimental values of the DWF 
for Ei = 63 meV and angles of incidence 19° (squares), 39° (circles), 51° (triangles), 
61.7° (inverted triangles) and 71° (diamonds). 

Fig. 4. Same as in Fig. 3 for Ei—21 meV and incident angles 31.8° (squares), 55.5° 
(circles) and 73.5° (triangles). 



Fig. 5. Values of the cut-off wavevector Q c as function of the normal component 
of He atom incident energy as determined from various procedures. Dashed line: 
Q c = (2a/z t f/ 2 where a is the range parameter of the Morse potential. This form 
of Q c has been used in the calculations leading to Figs. 3 and 4. Full squares: 
best fit Q c values for the vibrating repulsive component of the Morse potential. Full 
triangles: best fit Q c values for the adjusted exponential potential. Full circles denote 
the best fit Q c values for the vibrating Morse potential and they practically coincide 
with the dashed line for Ef > 15 meV. Open symbols denote Q c values pertinent to 
exponentially repulsive potentials given in the literature: open squares are from ||, 
open circle from ]I0| and open triangle from ||. 
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